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We propose a detailed description of the structure of the layer formed by polyelectrolyte chains 
adsorbed onto an oppositely charged surface in the semi-dilute regime. We combine the mean-field 
Poisson-Boltzmann-Edwards theory and the scaling functional theory to describe the variations 
of the monomer concentration, the electrostatic potential, and the local grafting density with the 
distance to the surface. For long polymers, we find that the effective charge of the decorated surface 
(surface plus adsorbed polyelectrolytes) can be much larger than the bare charge of the surface at 
low salt concentration, thus providing an experimental route to a "supercharging" type of effect. 

I. INTRODUCTION 

When a charged solid surface (charge density cr ) is exposed to a solution of polyelectrolytes carrying opposite ions (Z 
charges per chain), a counter-ion exchange occurs such that the chains adsorb and the counter-ions initially retained 
captive by electrostatic attraction are released in the solution. This is because the system gains the translational 
entropy of the micro-ions, which is ZksT (where ksT is the thermal energy), and only looses fc^T for each adsorbed 
chain. Then at sufficiently high charge density, adsorbed polyelectrolytes form a dense layer of coils strongly bounded 
to the solid. In some favorable cases, because the macro-ions have very long tails dandling in the solution, the adsorbed 
layer alone may carry a surface charge density larger than ctq. When the solution has been rinsed away, the decorated 
surface (solid surface plus attached chains) then behaves as a new charged system surrounded by counter-ions. This 
feature is the basic principle of a new device to build controlled charged multi-layers of alternated sign 1 . 

On the theoretical side, the structure of the adsorbed layer is known in the dilute regime, where the chains in the 
bulk are isolated coils, and the solution does not affect much the layer [2j |3j HI [5j |6j [3 [8] . Comparatively, much 
less is known about adsorption from semi-dilute solutions. Because the dilute regime for polyelectrolytes is found 
at vanishing monomer concentrations (4>* h ~ N~ 2 where N is the index of polymerization), the semi-dilute regime is 
certainly relevant from an experimental point of view. Note that the Decher process takes place in the semi-dilute 
regime in Ref . |9J . Borukhov et al. have numerically solved the non-linear mean-field system of equations in the semi- 
dilute regime for repulsive surfaces [10 . They also used scaling arguments to determine the relevant characteristic 
length scale in terms of the total electrostatic potential drop, but no detailed analytical description of the layer is 
done. Besides, the case of attractive surfaces does not simply follow from the repulsive case explicitly treated by 
Borukhov et al. since the concentration profiles arc significantly affected by the boundary conditions. Chatellier and 
Joanny have solved the linearized mean-field system of equations suitable to describe how the charged surface affects 
the semi-dilute solution [TT]. However, this approach is restricted to small perturbations. By solving numerically 
the self-consistent field equations using the ground state dominance approximation, Wang explored the influence of 
various parameters (surface charge, salt concentration, bulk concentration) on the occurrence of charge inversion |12j . 
The latter was found to be strong at high salt concentration. 

In any case, these works do not provide the extension, H, of the adsorbed layer (which does not identify with the 
characteristic relaxation length of the concentration), and the amount of material, T, attached to the solid surface 
(which does not identifies with the excess material driven to the interface J(4> — 0b) dz, where <j> is the monomer 
concentration, as defined in Refs. [TUJ [TJl US])- There is a technical reason for that : the mean-field approach for 
adsorbed polyelectrolytes in semi-dilute conditions has to be supplemented by another approach to provide the full 
picture. This is because the variations of the polymeric concentration and the electrostatic potential does not carry 
enough information to completely describe the layer. This matters in practical situation since T is directly related 
to the charge density carried by the decorated surface after removal from the solution. Accordingly, this quantity 
deserves a special attention. Moreover, it has been seen experimentally that the extension and the amount of material 
of the adsorbed layer vary significantly with N [14] . This effect is not explained using the mean-field theory in the 
ground-state dominance approximation [TS] alone, which is used in Refs. [TUl HH |T2J [IS] . Indeed within this approach, 
the semi-dilute solution is described with <f> and the adsorbed chains cannot be distinguished from free chains in the 
semi-dilute layer. 

In this article, we propose a complete description of the adsorbed layer in equilibrium with the semi-dilute solution. 
In particular, we estimate T and H as a function of N, a and (j)b- We use a mean- field approach supplemented by 
a scaling type of description. The mean- field theory that we consider is the celebrated Poisson-Boltzmann-Edwards' 
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(PBE) description of the polyelectrolyte solution (see the reviews [131 US])- The scaling approach that we use is the 
Scaling Functional Theory (SFT) [17] that we adapt to treat semi-dilute solutions of polyelectrolytes. 

We shall consider a rather standard situation: long (N ^ 1), linear, fully flexible (the persistence length is identified 
with the monomer size a), polymer chains, small fraction of charges (/ <SC 1), in 9 solvent conditions (t> = 0, and 
the monomer/monomer interaction is dominated by three body excluded volume interactions), bare surface charge 
(To. We assume no added salt and no specific interaction with the surface, as a first hint into the full problem. For 
definiteness, we suppose (without loss of generality) that the surface is negatively charged with a bare surface charge 
density —<j$ (with ctq > 0), and the chain carries Z = fN quenched elementary positive charges +e where / is the 
fraction of charged monomers. 

In view of the number of parameters, this issue is a formidable task. To proceed further, our analysis is based on 
two major assumptions : 

1. That we may find long polyelectrolyte tails belonging to adsorbed chains far from the surface (protruding in the 
solution), will be our first assumption (see Fig. |TJ) . This hypothesis is inspired by the results found with long neutral 
chains adsorbed from semi-dilute solution, where we know that the concentration field decreases over a distance 
comparable to the bulk blob size £b ~ <^ 3 ^ 4 (the characteristic length scale of the solution), whereas the layer of 
adsorbed chains extend over much larger distances H = Rq ~ N 1 / 2 ^ 8 [TTl [T5J [T2] . In this example, the essential 
idea is that the collective response of the chains screens the perturbation induced by the surface at a characteristic 
length scale which does not depend on the chain size [20]. A similar feature is found with polyelectrolytes: in the 
semi-dilute regime, the Debye-Hiickel length of the solution, kT 1 , is independent of the chain length. One may possibly 
argue that polyelectrolytes lie quite flat on the surface, as with dilute solutions. This is however very unlikely since 
for such picture extrapolated to very high concentrations of polyelectrolyte, it would mean that the adsorbed chains 
adopt an almost two-dimensional configuration. Clearly, the system of adsorbed chains would gain large amounts of 
conformational entropy by allowing more chains to adsorb and recovering an almost Gaussian structure (thus leaving 
some tails dandling in the solution), without perturbing much the distribution of charges in the vicinity of the solid. 
In other words, since kT 1 <C Rg ~ N 1 / 2 in the limit N ^> 1, we expect that the influence of the surface will be small 
for the chain structure in equilibrium. 

2. Looking at the variations of the electrostatic potential, ^f(z) with the distance to the surface, z, in the dilute regime, 
we anticipate a very simple form for ^(z) in the semi-dilute regime (see Fig. [2). Fr om the mean-field description of 
the layer in the dilute regime [2] [3], we learn that the electrostatic potential first increases from negative values, then 
saturates, and eventually relaxes to zero. That there is a regime of parameters in the semi-dilute regime where such 
profile for 'J will persist to some extent be our second assumption. This assumption will be justified by the calculation 
but one can find a simple argument for it. The monomer concentration profile in the dilute regime tells us that the 
osmotic pressure decreases from the solid surface. Suppose that we increase the concentration of the solution from the 
dilute regime into the semi-dilute regime. Because the osmotic pressure close to the solid surface is much higher than 
what is found in the solution, we expect that the presence of the solution will affect the structure of the outer fringe 
of the layer, not the structure of the inner regions. Accordingly, we do not expect that the form of the electrostatic 
potential will be profoundly affected, at least in the semi-dilute regime close to the dilute regime. 

Following our second hypothesis, we define the plane of vanishing electrostatic field (such that <N> /dz = 0) as the 
neutral plane. Similarly, we set the plane where the electrostatic potential vanishes (such that ^ = 0) as the relaxation 
plane (we implicit assume that we haven chosen ^(oo) = 0). Then the layer might be divided into three main regions 
(labeled (1), (2) and (3), see Fig. [2} such that: 

1. The compensation region stands between the solid surface and the neutral plane. By definition, the charges 
provided by the polyelectrolyte material situated in this region exactly compensate the surface charge. 

2. The intermediate region stands between the neutral plane and the relaxation plane. 

3. The outer region stands between the relaxation plane and the outer border of the layer, as defined by the longest 
polymeric coil in direct contact with the surface. 

It is important to realize that this distinction, although a priori, is not arbitrary. Each of these region is well defined 
from a physical point of view. 

In what follows, we first introduce the Scaling Functional Theory suitable to describe polyelectrolyte semi-dilute 
layers (Section 2). Then we propose a detailed picture for the adsorbed layer in Section 3. Section 4 presents 
our concluding remarks focusing on the issue of charge inversion. Scaling laws in polyelectrolyte solutions and the 
mean-field approach are briefly introduced in Appendices 1 and 2 respectively. 
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FIG. 1: Sketch of the adsorbed polyelectrolyte layer exposed to a semi-dilute solution. Some long dandling tails and loops 
are protruding in the bulk. The total monomer concentration profile is thus the sum of the adsorbed monomer concentration 
and the free monomer concentration. We note that the relaxation length of the monomer concentration close to the surface is 
different from the total adsorbed layer thickness, H. 
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FIG. 2: Sketch of the variations of the mean electrostatic potential in the direction perpendicular to the surface. It starts with 
a negative value at the negatively charged surface, increases until it reaches a maximum at the "neutral plane", where the 
electric field vanishes. Then it relaxes until it reaches its average value in the bulk, 0, at the relaxation plane and is constant 
in the outer region of the adsorbed layer. 



II. SCALING FUNCTIONAL THEORY 



This theory was first introduced to describe neutral polydisperse polymeric layers by generalizing the Alcxandcr- 
de Gennes model for brushes [21] . The basic principle is to describe the layer of adsorbed chains as a thermodynamical 
ensemble of tails and loops where the entropy associated with the polydispersity in size compete with the elastic energy 
and the monomer excluded volume. From a formal point of view, the SFT approach is a functional theory where the 
Hamiltonian is simplified according as follows [22]: a) loops are cut into pseudo-tails, b) the tails and pseudo-tails 
are described by a single function z, where z(n) is the position of the nth monomer, n, of each tail. This idea proved 
successful in describing various situation ranging from polymer brushes, reversible adsorption, irreversible adsorption, 
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whatever the solvent conditions. 

Our aim in this article is not to generalize the SFT to describe polyelectrolyte chains. Rather, we take advantage 
of the fact that in the semi-dilute regime, the chain is both Gaussian and neutral at large scales, and the system of 
polymers is thus analogous to a neutral solution provided that the blob is renormalized according to the physics at 
small scales. Essentially, our approach consists in integrating the degrees of freedom of the various species (including 
electrostatics interactions) into a proper description of the blob. 

Within the SFT, the system of loops and tails attached to the solid surface is described by three functions: S, 
the local grafting density of tails and pseudo-tails (a function of the curvilinear distance n), z, the path, and 4>, the 
volume fraction of monomers. These functions are linked through the conservation of monomers: 



S(n(z)) 
z(n(z)) 



For polyelectrolyte semi-dilute layers, the effective free-energy (per cm 2 ) of the system write ((3 1 = ksT): 




+ a-^/^Wi 3 / 2 ^) - S'(n) In 



So , 



(2.1) 



(2.2) 



where Sq is the "grafting density" at the surface and £ e ~ a(f 2 @B) 1 / 3 is the electrostatic blob size (in solvent 



conditions, £b is the Bjerrum length in units of a, cf. Appendix 1). The effective free-energy eq. ( |2.2[ ) is the sum of 
three contributions which account respectively for : 

a) monomer excluded volume interactions. This contribution is J Hdz, where the polymeric contribution to osmotic 
pressure, II, scales as ksT/£ 3 : 
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With eq. (2.1 ), it is simple to show that 
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where H — z(N) is the position of the last monomer and is identified as the layer thickness. 
b ) the elasticity of the tail at a scale larger than the blob size (after integration by parts) 
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c) the entropy associated to the polydispersity in size of loops and tails which is similar for charged and neutral 
monomers as soon as polymers are flexible |21j . 

Minimizing the free-energy eq. (2.2) with respect to S and z yields two equations 

1/2 



3<r-3/2 



s 



6 



3/2 



s 



3 

3/2 



£ 1/2 S 1/2 i 3/2 



3 

2n 



al /2 s 3 /-z 



2-1/2 



= n h 



(2.6) 



(2.7) 



These equations are formally similar to what we find in the neutral case (and their interpretation in physical terms is 
thus identical), except that the small scale structure introduces different scaling relationship. Equation (2.7 1 describes 
the local balance of forces (per unit surface parallel to the solid) : elasticity (first term of lhs) competes with osmotic 
pressure (second term of lhs) and the bulk pressure (rhs) . Note that the osmotic pressure of counter-ions is the same 
in the outer region and in the bulk so that it does not appear in eq. (2.7 1. Equation (2.6) describes the conservation 

fib = dUb I d(j>b is the chemical potential in the bulk and scales 



of the generalized chemical potential. In equation (2.6 1 
as or 1 : 



Mb - <P, 



1/2 



3/2 



(2.8) 



The bulk osmotic pressure, Ilf, ~ l/£ 3 , in equation \2.7\ accounts for the pressure induced by the solution. 



We emphasis that the SFT system of equations (|2.6p-(2.7l assumes that the properties of the tails and loops are 



that of neutral strings of blobs, and may only apply to situations where the electrostatics interactions are screened at 
distances larger than the mesh size £ and especially very close to the substrate. 
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III. STRUCTURE OF THE LAYER 

A. Compensation region (1). 

This region, close to the solid surface, is such that the charges provided by the polyelectrolyte material exactly 
compensate the surface charge. Hence 

Ti = * /f (3-9) 

and > in this region (at the outer boundary of the compensation region: 4^ = 0). For this reason, we expect that 
the counter-ions from polyelectrolytes are not present, p~ = 0. Moreover, at the low salt limit, positive counter-ions 
(from the charged surface) are negligible since their concentration is proportional to surface/ volume and vanishes in 
the thermodynamic limit. 

Each electrostatic blob carries fg e charges and occupies a surface £_ 2 . These relation defines a critical value for <tq 
such that a single "carpet" of close packing electrostatic blob is able to fully compensate the surface charge: 

o c = f^ 2 (3.10) 



The monomer volume fraction inside this region is thus </> e , given by Eq. (0.421. Moreover, the electrostatic energy 
of a blob interacting with the surface is of order fc^T for the critical value, oo ~ er c , i.e. of the same order of the 
electrostatic repulsion between monomers. Accordingly, we shall distinguish two cases: cto < a c and ao > r c (Fig. [3j) . 

III.A.l Case (T < Ve- 
in this regime, we rely on the results found by Dobrynin et al. [3] for the case of dilute solutions. The electrostatic 
attraction of charged monomers to the charged surface is weaker than the electrostatic repulsion between charged 
monomers. Therefore, the chain statistics is rod-like at scales larger than the blob £ (in the plane parallel to the 
surface) and Gaussian at smaller scales (£ ~ ag 1 / 2 ). However, the blob size is found by assuming that the electrostatic 

attraction to the surface is of order fc^T: ^ g<T °| e ~ ksT which leads to £ ~ £ e (fj) > £e- The chain segment is 

thus confined by the electrostatic attraction to the charged surface. However, contrary to the case of dilute solutions 
(single chain adsorption), polyelectrolytes do not lie flat on the surface since, as mentioned in the Introduction, the 
system gains conformational entropy by allowing more chains to adsorb with large loops protruding in the solution. 
The concentration of blobs in this region is found by assuming that this compensation layer actually compensates the 
surface charges and that the thickness is given by the confinement blob size £: 

/ \ 1/3 

Ai ^H^J (3 - u) 

This scaling result has been shown numerically for semi-dilute solutions by Borukhov et al. |10j and Shafir et al. |23] 
using the mean-field equations. The electric field vanishes at z — Ai in this regime. Of course, we do not have a 
closed packing of these blobs, £, since the number of blobs per unit surface is cro//<? — (oo/ctc) 5 / 3 ^ 2 < £e 2 - 

This picture is valid only if the average monomer volume fraction close to the surface is higher than the bulk value, 
i.e. e (cr o /CT c ) 4/3 > 4> b or 

3/4 

Ob ■ 



°c I ~J <<Jo<a c (3.12) 
III. A. 2 Case <t > a c . 

In the limit a er c , the monomer volume fraction close to the surface, <fio, increases being larger than <p e , the 
electrostatic attraction becomes larger than the electrostatic intra-repulsion, and the osmotic pressure tends to swell 
the comp ensation region. Assuming that excluded volume interactions dominate over elasticity, as in Ref. [3], solving 



eq. (0.48 1 of Appendix 2 yields (at the leading order in (j) e /(j)o <C 1 and assuming /i ~ g e 1 in this compensation region) 
a parabolic profile: <f>(z) = <j) e + (4>o — </> e )(l — z/Ai) 2 , where 4>o — (STriB&o) 1 ^ 3 anc ^ 



Al S £ e — (3.13) 



The structure of the compensation region is thus self-similar (see Fig. [3]) as already described for th e dilute case in 



Ref. [5]- For the general case, the polymer elasticity is no more negligible and the full equation (0.481 must be solved 
without any approximation. 
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FIG. 3: Sketch of the compensation region (1) in the 3 cases: oo < °"c 0"o — a c and o"o > &c- The thickness of this region, Ai 
decreases with ao for ao < a c and increases for ao > a c . 



It is important to note that our result Eq. (3.131 is in contradiction with Andelman et al. [TUl EH H3]. They 



predict that the characteristic length scale of the adsorbed layer decreases with ao ([13]) eq. (71)) following Eq. (3.11 1 
However, we argue that, for strong adsorption (cto S> <j c ), the structure of the compensation region is found by 
balancing the electrostatic adsorption with the osmotic pressure (excluded volume interactions) and that chains are 
not confined. Instead of blobs of confinement, the structure is described in terms of osmotic blobs, in a similar way 
of describing neutral polymer layers. 

Anyway, as we argue below (Section 4), in the semi-dilute regime and for long polyelectrolytes, we do not expect 
that this region has a dominant influence on observables such as the thickness of the layer, or the amount of material 
attached to the surface, and it is enough for our purpose to treat the restricted case where ao ~ a c , and 

Ai=£ e (3.14) 

Anticipating the use of the SFT, it is useful to estimate the amount of monomers (per tail), ni, involved in the 
compensation region. In a simple view where each tail contributes to one electrostatic blob, we find 

ni=g e (3.15) 

and the local grafting density at the outer limit of the compensation region scales as 

1 



S(rii) = „ ( 3 - 16 ) 



B. Intermediate region (2). 



In this region, the direct electrostatic influence of the surface vanishes, and from this point of view, the structure 
should be similar to a layer of polyelectrolytes chains "adsorbed" onto a Active neutral surface positioned at z = A\. 
Here, the concentration of monomers decreases from <p e to the bulk value, </V Because monomers and counter-ions are 
ruled by very different equations, the electrostatic potential ^(z) does not vanish. As the system is globally neutral 
and <^(Ai) — 4> e > (j>b in the semi-dilute regime, charge conservation in the system implies that the volume fraction 
of counter-ions, p~(z), will decrease from p~(Ai) > f(f>b to the bulk value f(f>b at the external border of region (2). 
Knowing that p~ is related to the adimensional electrostatic potential 'J = (3eV by = f<fib exp W (see Appendix 2) , 
we deduce that ^(Ai) = is positive (with the boundary condition \&(oo) = 0). Thus, we expect that ^(z) decreases 
monotonically from a finite positive value, ^ (to be estimated below), to 0. 

When (f>b is not too small in comparison to <p e , say <p e — cf>b <^ </>&, the linearized PBE system of equations described 
in Appendix 2 is justified to find the relaxation of the concentration profile. Note that for the case tfib <§C 4> e the full 
system should be solved numerically [TO], [23] . With the appropriate boundary conditions [<fi(Ai) — cf> e , <j)(oo) = <f>b, 
dV/dzikx) =0 and *(oo) = 0], we find: 



1-R 

4ir£ B f<l>b ~ a 2 A_ 



R 



-(*-Ai)/*4 



p(z) = /&[*(«) + !] 



(3.17) 

(3.18) 
(3.19) 
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where 



r ± 2 



\ = ^ (*l B f + H b ) [i ± VA] 

12^/(202 + /) 



R- L = 1 



A_ (4^ B / fc ~ a 2 A+) 



(3.20) 
(3.21) 

(3.22) 



In Fig. Hare plotted 5ip(z) — ip(z) — -/0b and ^(z) for £g = 1, / = 0.1, and 06 = 0.26 (concentrated semi-dilute 
solution)"! Note that the characteristic relaxation lengths £± were obtained by Chatellier and Joanny [TT]. In this 
reference, the possibility of damped oscillations in the volume fraction profile (£± imaginary) was carefully examined 
for lower values of 06 . Here, we will only consider the situation where £± are real numbers, for simplicity (the same 
discussion can be carried on with the real parts of £± in more general cases) . Hence the surface induce an external 
perturbation to tp of magnitude yftpl — v / 06 which relaxes, at the linear response level, to the bulk value, \/06 on the 
characteristic wave length of the fluctuations in the bulk. 




FIG. 4: Polymer order parameter (left), 8ip/yf4> b , and electrostatic potential (right), \f , /\£ r i, as a function of the distance to 
z = Ai (in units of £-), in the intermediate region for Ib = 1, / = 0.1, and (fib = 0.26, in the mean-field approximation. 



Of course, we identify A 2 with the longest relaxation length, £_ 

A 2 = £- 



(3.23) 



It is instructive to com pare t he characteristic relaxation length £_ with £&, which requires that we simplify the 
expression given by eq. (3.201. Assuming that the intermediate region is dominated by the counter-ions pressure 



with an effective excluded volume v c s ^ f /(fib (and that contribution arising from the monomers is negligible), the 
relaxation length £- writes (with A ~ 1): 



£- 



(3.24) 



In this limit, we have £_ < £b for 06 < (f /£b) 1 ^ 3 , & reasonable case when (fib < <fi e . The amount of monomers is this 
region is given by 



Mi+A 2 g 

- 3 / (fi(z)dz (0.1706 + O.430 e + 0.40 a/^^) 
Ja 1 a 6 



Equation (3.18) also yields an estimate for the extremum value of the electrostatic potential "Pi = ^(Ai): 

1 - R 



#1 = 8tt£ b 



/-v/06 ( V0e - V^bj 



R 



4tt^/06 - a 2 A_ (47rf B /06 - a 2 A+) 



Equation (3.26) in turn allows to estimate ^0 = = 0) through the integrated Poisson equation: 



*o = *i -47r^ B cr Ai +4ir£ B f / / (fi(z)dzdu 

Jo 



(3.25) 



(3.26) 



(3.27) 
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For the moderate adsorption case (oo ~ er c ) that we consider in this article, we find: 

f9e 



$0 = #1 ~ 

For the case of strong adsorption (<tq 3> <r c ), we obtain: 



2tt£t 



f9e ( &o\ 



4/3 



£e \<Jc) 



(3.28) 



(3.29) 



The SFT approach described in Section 2 is unlikely to describe correctly the intermediate region. At a simple scaling 
level, we will note that the local grafting density should drop from S(ri\) = A to S(ri2) = j?, where n 2 is the 
curvilinear index at z = Ai + A 2 . A crude estimate for ri 2 is found with a renormalized linear string approximation 
for the tails in the intermediate regime: 



A 2 



n 2 = ni+ g e 



9b 



which is equivalent to write 



A, 



(3.30) 



(3.31) 



where the last approximation holds in the limit (p e <C 4>b- Hence, the thickness of region (2) is then given by 
equation (3.231 at a mean-field level and for <f> e — <pb <C <pb, and by equation (3.31) at a scaling level. 



C. Outer region (3). 

Here, the polymer concentration (resp. the counterion concentration, the adimcnsional electrostatic potential) is 
the bulk concentration, <f>t, (resp. at scales larger than the bulk blob size), but the tails and loops in this 

region belong to chains in direct contact with the solid surface. At a scaling level of description, the layer is a close 
packing of semi-dilute blobs (size £). By analogy with the Rubinstein description of the bulk solution [53], we expect 
that electrostatic interactions are screened at distances larger than the mesh size £, and the system is amenable to a 
description in terms of neutral chains. Hence the SFT approach put forward in Section 2 is now appropriate. 

Obviously, free and adsorbed chains interpenetrate, and we expect a slowly vanishing concentration profile for 
the monomers which belong to chains in direct contact with the surface, as depicted in Fig. [5^. In what follows, 
however, we assume that the free chains do not penetrate in the layer of tails and loops in direct contact with the 
surface. Accordingly, we consider a sharp, Heaviside type of profile at the outer fringe (see Fig. ^p). A similar type 
of simplification is found in various theories of polymeric layers. In the context of polymer brushes in good solvent, 
e. g., Alexander and de Gennes assume an Heaviside step for the concentration profile in order to find the extension 
of the layer |25[ 126] . Their scaling result was confirmed by more refined theories where this constraint is relaxed 
(instead a parabolic profile is found) 27J . Similarly, Guiselin assumes no interpenctration in the context of molten 
solutions of neutral polymers adsorbed at an interface [2S]. Here again, the scaling result was experimentally confirmed 
by measuring the amount of material irreversibly adsorbed |29] , Because the semi-dilute solution is amenable to a 
description in terms of neutral blobs of strings of electrostatic coils [24, 30J, it is likely that the situation is analogous. 
In any case, we are confident that the scaling results that we find will not be dramatically affected by our simplification. 

The SFT equations suitable to describe this region are ( |2.6[ )-( [2~7| . In equilibrium, we expect that the bulk osmotic 
pressure exerted at the extremity of the layer [z = H) is balanced by the inner osmotic pressure, and the elasticity 



can be neglected in eq. (2.7). This yields 

S 



z 



(3.32) 



and thus [eq. (2.1 1], 4>{z) = 4>b- When this relation is introduced into eq. (2.6), this equation simplifies into 



t 1 / 2 on 

%S 2 + = (3.33) 



which admits a power law type of solution in the limit N — > oo, 

,/3/4 



Sin) * ai(Ua)^ (3 ' 34) 
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FIG. 5: We expect a slowly vanishing concentration profile for the monomers which belong to chains in direct contact with the 
surface (a). In our approach, we simplify this picture into a sharp, Heaviside type of profile at the outer fringe (b). 



Note that we find S(n 2 ) = £ b 2 ■ Now solving Eqs. (3.32 1— (3.34 1 (with the boundary condition z(n 2 ) = Ai + A 2 ) for 
z, gives a Gaussian structure at large scale: 



z(n) = Ai + A 2 + ai^/a)- 1 '^- 11 ^ 1 ' 2 



which readily yields [H = z(N)]: 



A-, = a 



N 1 ' 2 



(Wa) 1/4 <^ /4 



(3.35) 



(3.36) 



and 



r, = 



N ^ /4 A^/2 

S ( n ) dn - 2 (C I U/4 



(3.37) 



Note that the result eqs. (3.36 1-( 3.37 1 are precisely what we would find with scaling arguments assuming that the 
chains are Gaussian at large scale, and the SFT is only a formal way to recover them in this context. Figure [6] shows 
a scaling representation of the whole adsorbed layer in the case en — c c . On the other hand, 5* is directly related to 
the loop size distribution P through [22 



S(n) 



N 



P{u) du 



(3.38) 



and eq. (3.341 may be experimentally checked with a surface force apparatus [3~T] . 



IV. DISCUSSION 



The structure that we propose for the adsorbed layer clearly separates the influence of the charged surface which 
occurs only in the compensation region, from that of the bulk solution which act on the outer region, with an 
intermediate regime in between. We briefly describe the influence of 4>b an d &o on the layer structure and then discuss 
the occurrence of charge inversion for layers of very long charged chains. 

We expect that as we increase <pb from (f>l , we progressively build the outer region at the expense of the intermediate 
region (A 2 decreases), but without much consequences for the compensation region. This is because the osmotic 
pressure of the counter-ions (the dominant contribution to the physics in the intermediate region) is much higher near 
z = Ai than at z = A 2 . Arguably, any increase of the bulk concentration of counter-ions will predominantly affect 
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concentration electrostatic blob 




FIG. 6: Scaling representation of the adsorbed layer of polyelectrolytes in the case <jo — cr c in terms of electrostatic blobs (size 
£ e ) and concentration blobs (size 



the low osmotic domain of the intermediate region, and hardly the structure near z = A%, where the pressure is high. 
Of course, when fa reaches cj> e , the intermediate region disappears in our simple view, and the bulk solution should 
then affect directly the compensation region. In this regime where fluctuations are negligible, numerical mean-field 
approaches are well adapted |10| 112] . On the other hand, if we increase gq from er c , at constant bulk concentration, 
the surface potential |^ol increases and more and more monomers are involved in compensating the surface charge 
(n\ and Ai increase). Accordingly, the concentration profile in the compensating region is affected. But for very long 
chains (JV 3> U2), we do not expect that this should influence much the structure of the outer region. 

Note that our description of the compensation and intermediate region resembles, but is not quite similar to the 
description of the adsorbed layer from dilute solution proposed in Ref. [3] : the so-called self-similar adsorbed layer 
regime. Essentially, we keep the parabolic decrease of the monomer volume fraction close to the surface, until we 
reach <j) = <j) e (where neglecting the chain elasticity is no more justified), and then we use the mean-field equations 
(which account for the elastic contribution and the osmotic pressure of counter-ions) to find the relaxation of <fi below 
<j> e . In Ref. [3J, the region where 4> < <fi e is described in terms of electrostatic blobs. However such a description 
does not take into account the presence of the counter-ions in this region which is imposed by electroneutrality (the 
electric potential is positive). In technical terms, we make advantage of the fact that <f> relaxes to a finite value, (f>b, 
to determine the monomelic profile below <j) e by linearizing the mean-field equations. Because <$> goes to in dilute 
solutions, such calculation is not possible and would necessitate a numerical resolution. The difference between the 
two descriptions is not just a matter of refinement since the relaxation of <fi yields the variations of the electrostatic 
potential. With our approach, we provide both a coherent and complete picture for the variations of this quantity 



[eq. (3.18)] 



A crucial aspect of adsorbed polyelectrolyte layers is their ability to capture a very large number of counter-ions per 
unit surface, much larger indeed than . In particular, the driving force for the adsorption of oppositely charged 
chains depends precisely on the amount of potentially free ions. With this in mind, it is interesting to discuss the 
following situation. 

Suppose that we build an adsorbed layer from a semi-dilute solution, as described before. Then we remove the 
polymeric solution after the layer has equilibrated, but we suppose that the adsorbed chains remain attached to the 
charged surface, in a fashion similar to the scenario discussed by Guiselin for irreversible adsorption of molten liquids 
of neutral chains [25J. Finally, we put the decorated layer into a solution of pure solvent. The apparent charge (also 
called "net charge" in Ref. [3J of the layer in the dilute solution will be zero, because of the counter-ions present 
within (or in the proximity) of the layer. On the other hand, the decorated layer carries a huge amount per unit 
surface of potentially free counter-ions. Provided that we disregard the fluffy structure of the layer, the decorated 
surface effectively behaves as a new charged surface of " nominal" charge density 

Aa = a-a = f(T 1 + T 2 + T 3 ) - a 

= /r 2 + 4 /12 / 7/e >fe /4 ^ 1/2 (4.39) 

Hence, for long chains, we expect that Act may be much larger than ctq. 

Let us compare: a) a layer of polyelectrolytes reversibly adsorbed from dilute solution in equilibrium, and b) a layer 
of polyelectrolytes irreversibly adsorbed from semi-dilute solution put into a solution of pure solvent (the solid surfaces 
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have the same cro in both cases). To this aim, we need a precise vocabulary. We define three regimes, depending on 
the comparison of a with <7o and 2gq\ 

• a < cro: undercharging (cr = cro corresponding to charge compensation); 

• (To < <7 < 2<7o: overcharging (cr = 2cr corresponding to charge inversion); 

• a > 2oo: supercharging. 

For case a) the mean-field theory predicts exact charge compensation in the absence of added salt [2]. With added 
salt, Joanny finds overcharging, with possibly charge inversion in the limiting case of very high ionic force. In any 
case, there is no supercharging. It also appears that the surface controls the amount of adsorbed material in case 
a) (together with the ionic strength of the solution in case of added salt), not the chain length, nor the polymer 
concentration. A more recent analysis is provided by the work of Dobrynin et al. [5] who carefully combine the 
mean-field theory with scaling arguments. But this study does not change the conclusions concerning the regime of 
compensation. 

In contrast, for case b) our analysis strongly suggests that we have supercharging for sufficiently long chains, 
N > (cro//-r 2 ) 2 /[(£ s / 2 ) 1 / 6 0s /2 ], even in the absence of salt. 



V. CONCLUDING REMARKS 



In this paper, we describe completely the structure of the adsorbed layer from semi-dilute solution by identifying 
the physical mechanisms in each sublayers : electrostatic attraction by the surface charge vs. three body interactions 
in the compensation region (1), importance of the polymer elasticity and the counter-ions osmotic pressure in the 
intermediate region (2), and presence of long tails and loops of adsorbed chains in the distal region (3) at the bulk 
concentration. We give explicit expressions for the coverage including all monomers belonging to adsorbed polymers, 
r and the layer thickness H as a function of parameters, bare surface charge, bulk concentration which could be 
checked experimentally. 

Clearly, our description of the intermediate region suffers from several weaknesses. First, we consider a mean-field 
approach extrapolated outside its domain of validity. Secondly, we use a linearized theory, which is only justified 
for bulk volume fractions close to </> e . However, it is very likely that our essential scaling type of conclusions will be 
robust to any improved description of the intermediate region. 

Similarly, the influence of the solvent quality (unless it is bad solvent) is rather limited. This is because the 
electrostatic interactions dominates over the excluded volume interactions at large scales, and in actual facts, the 
solvent only matters in the compensation region and in the electrostatic blob. For example, it is simple to show that 
r 3 = a- 2 <f)l ,A '{lBf 2 ) 1/ti N 1 /' 2 and A 3 S a^ 1/4 {£ B f 2 ) 1/U N 1 / 2 in good solvent conditions. 

This study is limited to the very low salt limit. When the ionic force is very large, electrostatic interactions are 
screened at length scale larger than the Debye length and lead to an effective excluded volume v e i = f 2 /(2n s ) where 
n s is the salt concentration [2j[13]. The structure of the compensation is then modified but far from the surface, the 
theory of neutral polymer layer [17] can be applied by introducing the total excluded volume v + v e i (where v is the 
bare one). One find for instance T 3 = a- 2 (t> 7 b /S [v + f 2 /(2n s )] 1 / s N 1 / 2 . 

As far as we know, very few experimental measurements of the amount of material attached to the surface as a 
function of N have been done under the conditions we are interested (oppositely charged surface, weakly charged 
chains, no salt added, semi-dilute solutions). The experiments of Chibowski et al. on adsorption of polyacrylic 
acid (PAA) and polyacrilamide (PAM) to a Fe2C>3 solid surface, show a significant enhancement of both the layer 
thickness and the amount of material when the chain molecular weight is increased [14] . For instance, the thickness 
increases from 4.53 nm at molecular weight M w =170 kg/mol to 6.69 nm at M w =240 kg/mol for PAA at pH 
3 and for a polymer concentration of 500 ppm (semi-dilute regime). This increase is in a qualitative agreement 
with Eq. (3.361 and is related to the formation of larger loops and tails. Of course, the authors mention the 
important role of the pH, since dissociation equilibriums are monitored by pH and thus both the number of charges 
of the chains and of the substrate can vary. Such effect has to be included in the theory to do a quantitative comparison. 



Concerning the first step of the layer-by-layer process PQ, it should be pointed out that the adsorbed layer is 
globally neutral (d^ / dz(H) = 0) because counter-ions which penetrate into the layer compensate the decorated 
surface charge. Hence it is not the apparent charge which is important but the amount of free counter-ions which 
can be released during the adsorption of oppositely charged polyelectrolytes. This is a key point since, by assuming 
that adsorbed chains stay upon rinsing, the driving force of adsorption is the exchange of counter-ions related to 
the amount of potentially free counter-ions. The good observable is thus the amount of adsorbed polymers T, 
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and not the polymer surface excess relative to the bulk concentration. Hence the structure of the adsorbed layer, 
described in terms of loops and tails, is necessary to compute the coverage. We thus go beyond standard mean-field 
approaches [H |5J El [HI EH E2] which cannot distinguish between adsorbed chains and free chains in the distal region. 
As a result, we find that the issue of overcharging-supcrcharging relics more on the chain length N and bulk volume 
fraction <f>f, than salt concentration : charge inversion occurs even at low salt concentration which underlines the role 
played by the chain length as a tunable experimental parameter. Presumably, this feature opens the way to a both 
simple and efficient way to tune the size of the successive layers in the Decher process. Obviously, a closer look at 
the case where two layers of oppositely charged polymers interpenetrate and rearrange |32] is required before we can 
drive any conclusion on the Decher process |33| . 

We are grateful to L. Bocquet, E. Trizac, G. Decher and R.R. Netz for interesting discussions. 



APPENDIX 1 : SCALING LAWS 



In this Appendix, we briefly recall the basic scaling concepts for polyelectrolytes. We also provide the essential 
equations needed to understand the scaling arguments presented in the article. Our present understanding of weakly 
charged polyelectrolyte solutions without added salt is based on the concept of electrostatic blob [TBI ES] : because 
the different charges of the backbone repel, the electrostatic interactions will deform the otherwise isotropic coil into 
an elongated string. At a scaling level, the chain is pictured as a linear string of electrostatic blobs (size £ e , number 

of monomers g e ) such that /3 {f ^ =1. Let us assume O solvent conditions for simplicity (the two-body excluded 

1/2 

volume parameter v — 0). Then the Flory scaling law for random walks (£ e = agj ) yields 

£e = —TR and 9e = (0A0) 

(p£ B ) 1/3 (f 2 t B f 3 

a 2 

where £b — £^ is the Bjerrum length, in units of a. The extension of the polymer is 

L=-H e = aN(f 2 £ B ) 1/3 (0.41) 

and the local volume fraction (inside the electrostatic blob), 

e -(/ 2 ^) 1/3 (0.42) 

When the bulk concentration of monomers is larger than the overlapping concentration, ^ = NL~ 3 ~ A~ 2 , the widely 
accepted picture is that electrostatic interactions are screened at distances larger than the mesh size of the solution [24 . 
Accordingly, the solution is a close packing of semi-dilute blobs (size number of monomers gj,, fa = gfc such 
that the chain is a Gaussian walk of blobs at a scale larger than and a linear string of electrostatic blobs at a scale 
smaller than £&. We have [5U] 

,-1/2 , -1/2 

& = b n-e and 9b = b m (0.43) 

Eventually, when the bulk volume fraction exceeds <p e , the electrostatic blob vanishes, and the chain is Gaussian at 
all scales. 



APPENDIX 2 : MEAN-FIELD THEORY 



In this Appendix, we recall the equations of the mean-field theory used to describe polyelectrolytes (see, e.g. [TTl [T5] 1 . 
The basic idea is to describe the monomers as microscopic ions carrying a charge fe, with the following specific features: 
a) the translational entropy is negligible (this contribution scales as (<fi/N)log((j)/N) and vanishes for long polymers 
(N ^> 1)), b) an excluded volume interaction estimated in mean-field, c) an elastic contribution to account for the 
chain deformation. Then the energy (per cm 2 ) writes 



(3F 



1 

2a 



2 ^o*(0) 



- *p - 
[p+(lnp+-l)+p-(lnp--l)] 



a 
87r£i 
dz 



(W) 1 



dz 

^3 



(0.44) 



12 



13 



where the first term is the self-energy of the surface, the second term is the electrostatic energy of resp. the co-ions 
and the counter-ions and the self-energy of the electric field, the third term is the translational entropy of micro-ions 
(in volume fraction p ± ), and the fourth term accounts for the energy of the monomers (resp.: three body-interaction 



with amplitude w 



— „6 



Edwards elastic contribution, electrostatic energy of the charged monomers). Minimizing the 



energy yields the Poisson-Boltzmann- Edwards (PBE) system of equations: 



d 2 ^> 



dz 2 
p+{z) 
p~(z) 

a 2 1 d 2 y^(z) 
dz 2 



(z) = 



,-*(*) 



/^)+r W! -*W-(n t +M)e*W 



— n s e 
= (n. + M)e*« 

= <j>{z) 2 + f*{z)-n 



(0.45) 

(0.46) 
(0.47) 

(0.48) 



where n s is the volume fraction of salt in the bulk. Note that the mean-field theory is strictly valid for concentrations 
above <f> e for at least one reason: the average volume fraction, <fi, is not the local volume fracti on ex perienced by 
the monomers, <p e , in the semi-dilute regime, and the three-body interactions (the term <fi 3 in eq. (0.44)) is therefore 



incorrect. However, in the absence of any other theory to describe these systems, it is quite common to extrapolate 
the mean-field approach outside its strict domain of validity, and try to solve the PBE system of equation in the 
semi-dilute regime. Because \i is the chemical potential of the monomers in the solution, and this quantity is sensitive 
to the local concentration, it is tempting to correct for p = l/g e (<f> 2 in th eta s olvent), while keeping the mean-field 
estimate of the chemical potential (4>(z) 2 in our context) in the rhs of eq. ( |0.48 1. However, this approximation has a 
technical drawback since *&(z) = 0, 4>{z) = 4>b IS then not solution of the Edwards equation, as we expect far away 
from the surface. Here, we will keep the PBE system of equation with p, = (/) 2 . In Section 3.2, we use the Chatellier 



In terms of the polymer order parameter ip 



and Joanny approximation |llj where both the Boltzmann eq. (0.46 0.47) and Edwards eq. (0.481 are linearized (resp. 
around '3/ = 0, and <p = 
into 



the PBE system of equation simplifies 



(z) = 5— y/(j) b 8ip{z) + K H o~ W(*) 



dz 2 
p+(z) 

2 2 d 2 5<p 



(n s + /&) [! + *(«)] 



6\ 



dz 2 



(z) = A4>l /2 5 V {z) + f^{z) 



(0.49) 

(0.50) 
(0.51) 

(0.52) 



where <C 1, f — \/<h + Sip (with 6<fr <C v^b) an d « 2 = 9,Trt B n s /a 2 
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